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Abstract 

We consider a two-phase elastic solid subject to diffusion-induced phase 
transformation by an interstitial species provided by a reservoir. We de¬ 
rive a simple analytical model to quantify the effect of misfit strain on the 
kinetic of phase transformation and to calculate the amplitude of the well- 
know hysteresis cycle observed when a sequence of forward and reverse phase 
transformations takes place. 
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1. Introduction 

Metallic nanoparticles are characterized by fast hydrogenation and de¬ 
hydrogenation kinetics and hence are of particular interest for hydrogen- 
storage applications. A source of complication in the understanding of the 
adsorption/release kinetics in these devices is the stress accompanying hy¬ 
dride formation due to the misfit between the metal and the hydride lattice 
structures. 
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It is well known that the elastic misfits associated to phase transforma¬ 
tion give rise to a stress held that may affect considerably phase equilibria 
in multiphase elastic solids D3C3- In this respect, of particular importance 
is the analysis carried out by Schwarz and Kachaturyan [3], who examined a 
two-phase solid solvent in contact with a reservoir providing solute interstitial 
atoms. Their analysis shows that transformation-induced strain makes it im¬ 
possible for the two phases to coexist at equilibrium and that, moreover, it is 
responsible for the hysteresis loop observed in a cyclic adsorption-desorption 
processes; however, it does not address the issue of how misfit strain may 
affect phase-transformation kinetics. 

In order to investigate this issue, in this paper we address the kinetics 
of phase-transformation in a spherical particle in contact with a reservoir 
of interstitial atoms at prescribed chemical potential /j,r, as shown in Fig. 
1 below. General multi-held theories have been devised which describe the 



concomitant processes taking place in the solid, such as phase transformation, 
deformation, and diffusion mum- In order to arrive at an analytically- 
tractable model, we simplify the problem by assuming that there exists a 
sharp interface which separates the body in two phases: 

• a low-concentration a phase, occupying the region Q a and 

• a high-concentration /3 phase, which occupies the region f Ip. 

By changing the pressure (and hence the chemical potential) of the reservoir, 
adsorption/desorption of atom takes places, and, in order to accommodate 
the increased/decreased amount of diffusant, the phase boundary, assumed 
to be sharp and coherent, is set in motion. 
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We base our model on the continuum theory of elastic solids coupled 
with species diffusion developed in [7j. In particular, we use the notion 
of configurational force as a tool to describe the evolution of the material 
structure of a body — in this case, the phase interface. 

A crucial assumption is that chemical potential does not depart from the 
reference chemical potential fj, Q , where the two phases would be at equilibrium 
in the absence of misfit strain. This value of chemical potential represents the 
equilibrium state usually considered in thermodynamics, for which surface 
and elastic effects play no role. This fact provides us with motivation to 
assume that the Gibbs energy depends linearly on the difference between the 
local value of chemical potential and the reference one. A consequence of 
this fact is that concentration is prescribed on each phase. 

We find, as a result, that the evolution of the phase boundary is governed 
by a differential equation whose solution can be explicitly worked out, and 
provides us with a precise assessment of the characteristic time of adorption 
and desorption. Our framework is fairly general, and may be adapted to 
several applications of technological relevance. In particular, it may be used 
to characterize the kinetics of hydrogen adsorption in spherical nanoparticles. 

A similar problem has already been considered in [8], where hydride for¬ 
mation and dissolution was studied under the assumption of small departure 
of chemical potential from a reference value. The study in [8] is limited to 
nucleation. Here, we consider both phase nucleation and phase growth. Ef¬ 
fects of lattice strains on kinetics of hydration/dehydration have also been 
considered in [annum, under different modeling assumptions. 

It is also worth noticing that our model is intimately connected with that 
describing a coherent precipitate from a solid solution, a problem treated by 
Voorhees and coworkers in p2i 13] (see also the review |T], as well as the 
paper mi In their investigation of the role played by elastic effects on 
the growth kinetics, they consider that the precipitate grows in an infinity 
matrix. The boundary conditions for the diffusion equation in the matrix, 
which is an equation for the composition held, stipulated two values for the 
composition, a far held one and a interfacial one, the latter obtained from 
the Gibbs-Thompson condition. 

The aforementioned problem is amaneble to analytical treatment under 
certain simplifying assumptions, including that the precipitate is spherical 
and that diffusion is fast in comparison to the the interface motion. Thus, 
by considerind small departures about the incoherent equilibrium (no elastic 
and surface effects), a kinetic equation for the rate of growth of a spherical 
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inclusion of radius p in an infinite elastically isotropic matrix is obtained (see 
equations (75) to (77) of [Hj). This result can be obtained from the result of 
our paper by assuming that R goes to infinity, neglecting capilarity effects, 
and relating with the far-held composition. 

2. Setting the stage 

In this section we assemble the ingredients of the continuum theory of 
elastic solids undergoing diffusion and phase transformation expounded in 
[7], to which we refer those readers who may want additional details. 

We consider the three-dimensional elastic body pictured in Fig. 1. 
The body has the shape of a sphere of radius R , and is partitioned into two 
time-dependent regions: the region occupied by the a phase, and the 

region ftp(t), occupied by f3 phases. We denote by p(t) the radius of h2 Q (t) 
and by by ^(t) its boundary, namely, the evolving interface separating the 
two phases. We assume that such spherical symmetry be preserved during 
the evolution process. The fields that are relevant to our description are: 

• concentration c; 

• chemical potential /q 

• stress S; 

• displacement u; 

• linear strain e; 

• flux of diffusant h; 

• Helmoholtz free energy per unit volume -0; 

• Gibbs free energy per unit volume c o. 

Each field depends on the typical point x £ Q and on time t. Linear strain 
and displacement are related by the compatibility condition 0 

£ = symVu =: e(u). (1) 


1 Here sym denotes the symmetric part of a tensor. In components, (JT]) reads: = 

2 (uij + with a comma denoting partial derivative. 
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Stress and flux of diffusant must comply with the balance of standard forces 
and mass balance of diffusant, namely^] 

diver = 0, divh = 0. (2) 

We enforce the above-mentioned balance equations away from the interface 
5?. The corresponding balance statements at the interface take the form of 
jump conditions: 

H n = 0, [c]H- [h] -n = 0. (3) 

Here n is the unit vector field pointing in the radial direction; V is the velocity 
of the interface in the direction of n; double square brackets enclosing a field 
denote the jump of that field at the interface. The first of ([3]) encodes the 
requirement that the normal traction be continuous across the interface; the 
second, that there be no production/depletion of diffusant at the interface. 
In addition to (J3|, we enforce at the interface the following conditions: 

N = o, [/i] = 0. (4) 

The jump conditions ([ 3 ]) state that both displacement and chemical potential 
are continuous at the interface. Such requirements encode our expectation 
that the interface be coherent , that is to say, that the displacement field does 
not jump (and hence no crack appears) at the interface. 

Continuity of chemical potential at the interface is known as assump¬ 
tion of local chemical equilibrium. The conditions Q at the interface are 
accompanied by0 

ern = 0, (jl = hr (5) 

at the boundary of the specimen. 

The Gibbs free energy per unit volume is determined by specifying two 


2 Using Einstein’s convention for repeated indices, the balance equations © are ren¬ 
dered componentwise as follows: 


'<■ 0,3 


= 0 , 


hj,j 0 . 


3 Operating pressures for hydrogen-storage in metallic hydrides are of the order of 10 6 — 
10 7 Pa. Such pressures, by themselves, produce modest strains and can be neglected. 
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functions co a (/i,e) and Qp(n,e) such that the constitutive equation^ 
to = £)</,(//, e) in for 4> — a, /?, 


( 6 ) 


holds at all times. Once the constitutive mapping delivering Gibbs energy 
has been prescribed, concentration and stress in phase are given by, re¬ 
spectively, 


c = 


e) 


and 


(T = 


du^n, e) 


d/j de 

As to the flux of diffusant, we assume it to have the form: 


( 7 ) 


h = —m^V/i, (8) 

where denotes the mobility of the diffusant in phase (j) = a,/3. 

In the next two sections we shall consider two issues: the first, is the 
determination of the fields of interest given the position and the velocity of 
the interface; the second is the determination of the extra condition that 
provides us with the equation of motion for the interface. 


3. The Gibbs energy 

Let us denote by g(t) the distance of the interface from the center of the 
sphere, and by g(t) its derivative, both at time t. In principle, when initial 
condition for concentration are prescribed, the evolution of g determines 
uniquely all fields of interest. Yet, the problem appears too complicated to 
extract useful information through an analytical approach. 

On the other hand, forward and reverse transformation take place in 
a neighborhood of the transition chemical potential po at which the two 
phases would be in equilibrium in the absence of elastic misfit. This fact 
suggests that we restrict ouselves to situations where the chemical potential 
/j differs modestly from /i 0 . Accordingly, we may assume that the Gibbs 
energy depend linearly on fi — /i 0 and take: 

= -04,(11 - / j , q ) + W^(e) (9) 

where the strain energy W ,depends on the particular phase. An immediate 
consequence of (j9]), which is crucial to our developments, is that: 


4 Equation ([6]) is read as follows: given a point x £ Q and a time t such that x is not 
part of Y(t),: ui(x,t) = cu a (y(x, t), e(x, t)) if x £ and ui(x,t) = ujp(y(x,t),£(x,t)) 

if a: £ Qp (t). 
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concentration is 


• in the bulk — that is, away from the interface — 
constant. Precisely, we have0 

c(x,t) — Cfj, for i e O^(t). (10) 


As a result, the conservation of diffusant reads divh = 0 and hence, in view 
of ([8]), chemical potential is harmonic in the bulk , that is to say, it solves the 
Laplace equation: 

A/i = 0 (11) 

away from the interface. 

On taking as reference the state in which the material is in the a phase 
and on assuming isotropic linearly elastic response, we are led to the following 
choice for the strain energy^ 


W t (e) 


G\e\ 2 + ||tre| 2 if 0 = a 

G\e — e 0 j' J + ||tr(e — £ 0 )| J if <t> = P, 


( 12 ) 


Here, G and A are the standard Lame moduli, assumed for simplicity to be 
equal in either phase and 

• e 0 is the misfit strain between a and f3 p/iaseQ 

Consistent with the isotropy assumption on the strain energy, we take the 
misfit strain to be proportional to the identity tensor ifl 


£ o — e oI, 


(13) 


through a constant e 0 . 


5 Eq. (101 means that if x £ P Q (t) then c(x,t) = c a ; if, instead, x € P^(t) then 

C(z,f) = C/3- 

6 Here |e| denoes the norm of the strain tensor, while tr is the trace operator. 

'In general, the misfit strain would depend not only on phase, but also on concentration. 
In particular, a concentration gradient may give rise to a stress for a single-phase specimen. 
In the present case, however, concentration is prescribed in each phase because of (101. 


8 In Cartesian components, (131 reads (eo)ij = eoA 


ij- 
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4. The boundary-value problems governing displacement and chem¬ 
ical potential 

We now can summarize the boundary-value problems governing the pri¬ 
mary fields when the position g(t) of the interface and its velocity V(t) = g(t) 
are known. The boundary, value problem governing displacement u and 
stress a is: 

diver = 0 

cr = 2Ge(u) + A(divu)I 
cr = 2Ge(u) + (A(divu)(2G + 3A)e 0 )I, 

H = 0 

W n = o 

(Til = 0 

The boundary-value problem governing chemical potential p and flux of dif- 
fusant h is: 

divh = 0 

h = —ttTq, V /i 
h = —mpS7p 
M = /hn 

[h] ■ n = (c 0 - c a )V(t), 

g = fin 

On account of the spherical-symmetry features of domain and data, one can 
look for solutions of the form: 

u(x,t) = u(r,t)n(x), p(x, t) — v(r, t). (16) 

where r = r{x) = |x is the distance of point x from the center of the sphere 
and n(x) = Vr(x) is the radial vector field of unit norm. For the sake of 
brevity, we omit the analytical expressions of u and v. Indeed, as we shall 
see in the next section, the only information we actually need to determine 
the evolution of the interface consists of: 

• the value of chemical potential at the phase interface: 
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in £} a (t) U Qpft), 
in Vl a (t), 
in Qp(t), 
on y{t), 
on y(t), 
on <90. 


in 0 Q (f) U O^ (f), 
in O a (t), 
in O p(t), 
on y(t), 
on y(t), 
on (90. 


(V) 



• the jump of strain energy at the phase interface: 

W,(s) - (2G - 3A - SgQ 3 ) 4; (18) 

• the jump of the normal derivative of the displacement gradient at the 
interface: 


2 G + 3A 

|Vu|n = lGTT £ » n; 

the value of the radial stress at the interface: 


1 \ ^ + 3A 

= 3 C 2CTI 


-'I 


(19) 


( 20 ) 


To keep the present contribution reasonably short, we omit the details of 
the calculations leading to the above set of equations. These details shall be 
provided in [15J. 


5. Configurational balance 

We may summarize the discussion in the previous sections in the following 
statements: 


if the radius g(t) of the a phase, and the velocity V{t) = g(t) of the 
phase front are known, the boundary-value problems (14)- © can 
be solved explicitly to obtain the instantaneous values of the fields 
of interest in the bulk; 

an additional condition is needed in order to prescribe the motion of 
the interface. 


When seeking the aforementioned condition, it is important to bear in mind 
that the phase interface is composed of different material points at differ¬ 
ent times; therefore, it does not possess intrinsic material identity. Thus, 
the evolution of such material structure cannot be described by standard 
Newtonian forces, whose power expenditure is associated to the motion of 
individual material points. 

In order to construct models that can describe the evolution of material 
structures, the point of view adopted in [T5] has been proven to be expedi¬ 
ent. The main idea is to introduce forces that expend work on the motion 
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of these structures and are subject to their own balance equations. In mod¬ 
ern continuum mechanics, these forces are called configurational or material. 
This concept has been succesfully used to model the evolution of material 
structures such as phase interfaces, cracks, in three dimensional as well as 
lower dimensional continua, such as beams mt 

As extensively discussed in PTB], the law governing the motion of a sharp 
interface can be derived from a configurational balance [C]n + g y + e y = 0 
involving a configurational stress: 


C = ul- VuOr. 


( 21 ) 


with g y and e y , respectively, the internal and external configurational force. 




It is worth noticing that the configurational stress, as defined in (21), is 


intimately related to the tensorial quantity introduced by Eshelby in his 
investigations on the motion of defects in elastic solids [IS]. 

For a structureless and dissipationless interface, both the internal and 
the external configurational forces vanish, and the configurational balance 
reduces to: 

[C]n = 0. (22) 


On taking the dot product of both sides of (22) with n, and on accounting 


for the continuity of traction and chemical potential across the interface, we 
arrive at: 


- ic :> - Ca)(p - Aid) + Wp(e) - W a (e) - [Vu]n ■ Sn = 0. 


(23) 


On substituting (IT)—(20) into (23) we obtain one of the main result of this 


paper, namely, the evolution equation for g(t), the radius of the a phase: 


Q 1 \ C/3 Ca 

n ~ d ~^ ee 


Hr \ , 0,2^-^ + 3A 


= (*-c °)^^~TJ +2y ° G 2G + X 




(24) 


6. Transformation time 


With a view towards extracting information from (24), we introduce the 
characteristic time: 


R 2 


T = 


mp/i o 


(25) 
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and we introduce the dimensionless quantities: 


V = 


r: 


a = 


= b = 


2 G (2 a + 3A)e, 


A*0 


2G + A ^o( c /3 ~~ c a ) 


Then, (24) takes the following form: 

TJ/(1 - y)y = 1 - a +1>( 2 t/ 3 - 1 ). 


(26) 


(27) 


We now look for solutions of (26) such that £>(0) = 1 and g < 0. Any such 
solution describes a physical process in which the specimen is initially in the 
a phase, the (3 phase nucleates at the boundary, and the phase interface 
moves towards the interior of the particle. This process is accompanied by 
adsorption from the reservoir until the specimen is saturated, a situation that 
corresponds to g — 0, when the particle is entirely in the f3 phase. We refer 
to this process as the a —> j3 transformation. It turns out that: 

• the a —> /3 transformation can take place only if the following condition 
is satisfied: 

a > 1 + b] (28) 


moreover, if (28) holds, then the transformation time is: 

T a ^n = — A/( k ). 


(29) 


where the parameter k and the function / are defined by, respectively, 


k = 


and 


/(«) = 


*(1 


! 0 1 + k(2z 3 — 1 ) 


dz. 


In a similar fashion, we can consider a solution such that £>(0) = 0 and g > 0. 
This solution describes a physical process in which the specimen is initially 
in the the f3 phase, the a phase nucleates at the center, and the phase front 
moves towards the boundary until the specimen is completely in the a phase. 
It turns out that: 

• the (3 —y a transformation can take place only if a < 1 — 6; in this case, 
the transformation time is: 

= T ^/(«). (30) 

1 — a 
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/(«) 



Figure 2: Plot of f(n). Positive values of k yield the transition time for the /3 —> a 
transformation. Negative values of k yield the transformation time for the a —> /3 trans¬ 
formation. Remarkably, elastic misfit slows down the [3 —* a transformation and speeds 
up the a —> j3 transformation. 


It is worth noticing that / has a closed-form expression, as can be verified 
using standar computer algebra systems. The graph of / is plotted in Figure 
2. Positive values of n are relevant in the a —> (3 transformation. Negative 
values are for the [3 —» a transformation. From (29) and (30) we can draw 
the following important consequence: 


• the ratio between the transition time with and without misfit strain is 
given by /(«)//(0). 

On taking into account the above observation, and on looking at the plot 
of / for k E [0,1), we reckon that the elastic misfit between the two phases 
substantially slows down the transition from the [3 to the a phase, with the 
tranformation time blowing up as k —> 1. Concerning the a —» [3 transfor¬ 
mation, we observe, on the other hand, that an increase of the elastic misfit 
yields a decrease of k and in turn, as is apparent from Figure 2, a reduction 
of the transformation time. We therefore conclude that the presence of a 
lattice misfit accelerates the transformation from the a to the f3 phase and 
plays a crucial role in the kinetics. 

As a final remark, we observe that in order to trigger the a —> [3 transfor¬ 
mation and the (3 —> a transformation the chemical potential must satisfy, 
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respectively, fi R > /i 1 / ' /io(l + b) and /i R < := /i 0 (l — b). In par¬ 

ticular, if the chemical potential of the reservoir oscillates between and 
\x\ ^G a , then the path of y versus y R would define an hysteresis loop, consistent 
with the results in [3|. 

7. Conclusions 

We have considered the continuum theory proposed in [7j describing dif¬ 
fusion of a chemical species in an elastic solid partitioned in two phases 
separated by a sharp interface. We have specialized the theory to a spher¬ 
ical domain containing a structureless and dissipationless concentric phase 
interface, and we have investigated the effects of phase-trasformation strain 
on the kinetic of the interface. 

A key point of our approach is the assumption that the Gibbs energy de¬ 
pend linearly on the difference between the local value of chemical potential 
and a reference value. The main consequence of this assumption is that con¬ 
centration is constant in each phase. This state of matters makes the ensuing 
mathematical problem tractable, for it allows us to compute explicitly the 
instantaneous values of displacement and chemical potential in terms of the 
position of the interface. When this information is fed to the configurational 
balance equation, we obtain an autonomous first-order differential equation 
governing the position of the interface. Our analytical treatment of this 
equation shows that the presence of a lattice misfit between the two phases 
accelerates phase transformation processes which start with one phase nucle¬ 
ating at the boundary, and slows down processes that begin with nucleation 
of one phase from the center. Furthermore, our analysis confirms that elastic 
mistfits, besides being responsible for hysteresis in the phase-transformation 
process, strongly affect its kinetics. All these effects can be quantified in 
terms of a modest number of parameters. 
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